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i— h ■ Abstract 

In this paper we study the limiting distribution of the k smallest gaps be- 
tween eigenvalues of three kinds of random matrices - the Ginibre ensemble, 
the Wishart ensemble and the universal unitary ensemble. All of them fol- 
low a Poissonian ansatz. More precisely, for the Ginibre ensemble we have 
a global result in which the k-th smallest gap has typical length n 3 / 4 with 
^"J . density x 4k ~ 1 e~ x4 after normalization. For the Wishart and the universal 

| unitary ensemble, it has typical length n~ 4 / 3 and has density x 3k ~ 1 e~ x after 

■ normalization. 



> 1 Introduction 

o 

(NJ ■ The goal of this paper is to establish the limiting distribution of the k smallest gaps 

between eigenvalues. In particular, we focus on the following three ensembles: the Ginibre 
ensemble (where the entires are i.i.d. complex Gaussian), the Wishart ensemble (where 
it is the sample covariance matrix of n independent complex Gaussian vectors), and the 
universal unitary ensemble (where the matrix is invariant under unitary transformations 
with a general potential). More precisely, we address two problems here: pick an n x n 
matrix from a given ensemble, what is the typical length of the k smallest eigenvalue 
gaps? As a second order approximation, under the correct scaling, what is the joint 
limiting distribution of these gaps as we send n to infinity? 

Through many years typical spacing between eigenvalues of random matrices has 
attracted a lot of attention, this is partially due to Montgomery's conjecture: the nor- 
malized gaps between zeros from a Riemann zeta function appear to share the same 
statistics of the normalized eigenvalue gaps of a circular unitary ensemble (CUE) (see 
[12] for a review). We already know that in the Gaussian unitary ensemble (GUE) case, 
the expectation of the empirical distribution of the normalized gaps has the limiting den- 
sity which can be expressed as a Fredholm determinant. p(s) := d 2 ss det(I — K,)^^ ^ where 
/C is the integral operator on L 2 (0, s) with the kernel K(x,y) = sin[7r(x — y)]/[7r(x — y)]. 
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The asymptotics of p(s) when s — > and s — > oo have also been obtained (see Chapter 
6 and Appendix 13 of [B]). Later Kats and Sarnak in [T3] proved that this can also be 
extended to the CUE case and the corresponding orthogonal and symplectic case, with 
convergence rates provided. Similar results for universal unitary ensembles (UUE) are 
proved in [TJ] . For the second order asymptotics, Soshnikov established in [15] the Gaus- 
sian fluctuation around the limiting distribution for the CUE case. Recently, Tao and 
Vu proved the Four Moment Theorem in [16J saying that local eigenvalue statistics are 
universal among all Wigner matrices provided that the first four moments are matched. 
Yau et al also proved the similar universality result (under different assumptions) by 
analyzing the time-reversal approximation of the Dyson Brownian motion. See [17J and 

m. 

However, little was known about the extreme spacings, even in some simple ensem- 
bles. Vinson in [I] first studied this subject and obtained results of smallest spacings for 
the CUE and GUE case by asymptotic moment estimation. More precisely, by matching 
the moments, he proved that the number of eigenvalue gaps in (0, n -4 / 3 s) tends to the 
Poisson distribution. Later Soshnikov in [5] established similar results for general de- 
terminantal random point processes on the real line with a translation invariant kernel. 
Although the main idea is also moment (more precisely, cumulant) matching, he used 
another technique: he introduced the s-modified point process and investigated the cor- 
relation function of it instead of the original process. Recently Ben Arous and Bourgade 
in [1] developed Soshnikov's methods to solve the problem for CUE and GUE completely, 
showing that the joint distribution of the k smallest gaps has a Poissonian limit. For 
other ensembles, there are almost no results concerning these. In this paper, we extend 
Ben Arous and Bourgade's method to analyze the smallest eigenvalue gaps of the Ginibre 
ensemble, the Wishart ensemble and the general universal unitary ensemble. 

For the Ginibre ensemble, we have two major issues. First, in the CUE or GUE 
case considered in [I], the eigenvalues (or eigenangles) are real. Hence we have a natural 
definition for the gap, i.e. Aj+i — Aj, thanks to the linear ordering in R. In our case, 
however, the matrix is not hermitian so we have to deal with complex eigenvalues. Of 
course for any eigenvalue Aj it is natural to define min^ \Xj — A» j to be the gap. However 
the point process J^™ =1 <5 n 3 / 4 min j7 ^ cannot have a Poissionian limit because there are 

at least two identical points. The key problem here is that the search direction must be 
a half plane — you can not "look back" when searching for the gap (recall in the real 
case we define the gap to be A i+ i — Aj rather than min{Aj — Aj_i, Aj + i — A;}). Thus our 
remedy is to introduce a total ordering in C and then consider the modified point process 
n 3 / 4 iam j ^ i {\\ : j-x i \,\jy\i} instead. We prove that the modified process has a Poisson 
limit. See Theorem 11.11 

However, the total ordering appeared in Theorem 11.11 seems cumbersome, and our 
main concern is the joint distribution of the k smallest elements ti,...,tk of the set 
{\Xj — Aj| : % 7^ j} — with no ordering involved at all! We prove in Corollary II. II that for 
the small gaps the ordering introduced above will have hardly any effect. More rigorously, 
denote t\, . . . to be the k smallest elements of {\Xj — Aj| : Xj y Aj}, then we prove 
that (t\, . . . ,tfc) = (ti, . . . ,tk) with high probability, enabling us to analyze the gaps of 
the natural version from the "ordered version" . 

To state the second issue, we start from the GUE case. From the semi-circle law, 
the empirical measure of the eigenvalues has a semi-circle density which is supported 
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on [—2,2]. For any e > we refer (—2 + e, 2 — e) as the bulk of the eigenvalues and 
(—2 — e, — 2 + e) U (2 — e, 2 + e) as the edge. Ben Arous in [TJ stated the result only for 
gaps in the bulk. This is due to a technical issue of the Plancerel-Rotach asymptotics 
of Hermite polynomials. For our Ginibre case, things are similar. Due to the circular 
law, this time the limiting measure is a uniform distribution on the unit circle in C. The 
kernel function K(x,x) = ir~ 1 e~ n ^ Yji=o( n \ x \ 2 Y/® converges to 7r _1 if |x| < 1 and to 
if \x\ > 1. But the convergence rate is slow for x near the boundary. Nevertheless, our 
result for the convergence of the smallest gaps is global, that is, we do not restrict the gaps 
to be within the bulk of the unit circle — they are freed not only to the edge, but further 
to the whole complex plane. This is done by a more careful analysis of the convergence 
rate of the kernel function, as is stated in Lemma I3TT1 Meanwhile, in the following lemmas 
we analyzed separately the inner part, the boundary part and the outer part of the unit 
circle to eliminate the edge issues. In our situation, the orthogonal polynomial in the 
complex plane is just z n . Compared to the Hermite polynomial in the GUE case, this is 
much simpler, which enables us to perform a more precise analysis. 

It's generally accepted (for example see [6]) that the eigenvalue gaps are anti-correlated 
with each other, that is, large eigenvalue gaps tend to be followed with small ones and 
vice-versa. However, our result that the small eigenvalue gaps follows a Poissonian ansatz 
seems to imply that they are asymptotically independent. This has a natural explanation. 
As the gaps are anti-correlated, the small gaps are very unlikely to cluster together. That 
is, for the real case (such as GUE), eigenvalues in the set {Aj : Aj+i — Aj < const ■ n~^ 3 } 
must fall apart, the distance between these Aj's are large enough so that these small gaps 
Aj+i — Aj are nearly independent. 

The rest of the paper is organized as follows. In subsection 1.1, 1.2 and 1.3, we will 
state the problem and our result, respectively, for the three ensembles. In section 2, we 
will provide a general overview of our methods. The proofs for the three ensembles will 
be given in section 3 to section 5, respectively. In section 6, a conclusion of the paper 
will be given. 



1.1 The Complex Ginibre Ensemble 

Let X n = (xij)2j = i G C nxn be a random matrix with independent standard complex 
Gaussian entries. That is, 

XlJ = 7! " 

where independent real Gaussian random variables with mean zero and vari- 

ance one. Here we normalize each entry by the factor a/2 to keep the variance unity. Such 
X n is called Complex Ginibre Ensemble. For further discussion about this ensemble, see 

®- 

Define A n = X n j ^fn as the normalization X n . Furthermore define Ai, . . . , A n as the 
n eigenvalues of A n . By the circular law (see [7]), the empirical spectrum distribution 
will almost surely converge to the uniform distribution of the unit disk in the complex 
plane. 

Since we are dealing with complex eigenvalues now, we need to define a total ordering 
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in C first. For two complex numbers Z\ and z%, we define 

zi ^ z 2 w | = ^ < <±.±; 

It is obvious that (• -< •) is just the standard lexicographical ordering. 

Without losing generality, we can assume Ai H A2 di • • • di X n . Consider the point 
process on R + x C 

n— 1 

X^EW.-a.IA) (1.2) 

i=l 

where = argmin^ =1 • -^{|Aj — Aj| : Aj >z Aj} is the index such that Aj* is the closest one 
to Aj among all the eigenvalues greater than or equal to Aj. 

The main result is about the convergence of x to a Poisson point process. 

Theorem 1.1 As n — > 00, the process x converges weakly to a Poisson point process 
X in IR + x C with intensity 

E x{ AxI) = ([Wdu)([ ^V RnM />*. 

\Jb J V/nD(0,l) 71 J 71 J A 

for any bounded Borel sets A C M + and I C C. Here D(z,r) is a disk centered at z and 
having radius r. Moreover B = {u G C : \u\ G A, u >z 0} and | • | denotes the Lebesgue 
measure of the set. 

Remark 1.1 In the definition of i* = argmin^ =1 • -^{|Aj — A«| : Xj >z Xi}, the constraint 
Xj >z Xi is necessarily. Suppose A^ and Aj 2 are the two closest eigenvalues. If we didn't 
impose the constraint, we would simultaneously have i\ = %i,i\ — i\ and hence |Aj — Aj*| = 
\\.x — Xq\. This imply that at least two points in x share the same first coordinate. 
Thus x^ cannot have the Poissonian limit. 

From Theorem 11.11 the intensity is proportional to \I PI D(0, 1)|. This is intuitive 
because by the circular law, the limiting empirical distribution is uniform inside the unit 
disk and the region outside D(0, 1) does not contribute to the intensity. 

As a comparison to non-repulsive i.i.d. samples, we define fa, 1 < i < n to be 
independently and uniformly distributed in the unit circle. We can define a similar point 
process 



£(n) = ^<J (n | .._ .| i0 .). 



i=l 



where i* are similarly defined as above. We have the similar version of Theorem 11.11 

Theorem 1.2 As n — > 00, the process x^ converges weakly to a Poisson point process 
X with intensity 

m A x !) = ^M! / rd r. 

K J A 

for any bounded Borel sets A G M. + and JcC. 



4 



We omit the proof of Theorem 11.21 which is much simpler compared to Theorem 
11.11 We observe that the typical length of gaps between the eigenvalues from the Ginibre 
ensemble is 0(n~ 3 ^), while for i.i.d. case it is 0(n~ l ). The larger gaps in the Ginibre 
ensemble demonstrates the repulsive force. Even after normalization, the intensity cor- 
responding to Ginibre ensemble has the r 3 dr term, exhibiting heavier tails compared to 
rdr in the i.i.d. case. 

Now let's return to the main theme of the paper. Let < ... < t[ n) be the 
k smallest elements in the set {|Aj — \j\ : i 7^ j}. Moreover let rj^ = (7r/4) 1 / 4 ^ n ' ) , 
I = L2 . . . , k. We have the following corollary which describes the limiting distribution 

oiT k ■ 

Corollary 1.1 For any < X\ < y\ < . . . < x^ < yt, o,s n — > 00 

P(x e < n 3/4 r, (n) < y t , 1 < £ < k) — * Lr< - e~A \[{yt - x 4 ). 

^ ' i=\ 

In 'particular, n 3 ^ 4 r^ n converges in distribution to t\. with density 

P(r fc e dx) oc x^e'^dx. 

For a special case k = 1, it describes the distribution of the minimum gap between 
eigenvalues of the Ginibre ensemble. The density is proportional to x 3 e~ x . 



1.2 The Complex Wishart Ensemble 

As the second part of the paper, we consider the complex Wishart ensemble. Let the 
matrix X mn = (xij) G C mxn be a rectangular random matrix of Ginibre ensemble. That 
is, x^ = (uij + \J— lt>y)/v2 where Uij,Vi/s are i.i.d. standard normals. A complex 
Wishart ensemble, denoted as W 2 (m,n), can be obtained as A^ n A mn . Here X* denotes 
the conjugate transpose of X. Since X^ n A mn and X mn X^ n share the same non-zero 
eigenvalues, we can assume without losing generality that m > n. 

Define A mn = X^Xmn/m e c nxri as the normalization of the Wishart ensemble. 
Such form can be recognized as the sample covariance matrix. Denote Ai < . . . < A n 
as the eigenvalues of A mn . As we send n, m — > 00 while keeping the ratio fixed m/n = 
j5 G [1, 00), the empirical distribution of these eigenvalues converges in probability to the 
Marchenko-Pastur Law (see [S]). The density is 

f x l_ ^{{i + p-^y-x){x-{\-^W) 

9[X) ~ 2vr ' x 

Similar to the Ginibre case, for any pre- fixed eo > 0, we define the random point 
process in IR + x R by 

n-l 

x (n) = ^5 {n 4/3| Ai + 1 _ Ai | iAi) l {(1 _ /3 -l/2 ) 2 +eo<Ai<(1+( g-l/2 ) 2_ eo} 

i=l 

Note that here we do not need to define a total ordering because there is already a natural 
one on R. 

Our second result is stated in the following. 
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Theorem 1.3 For any fixed eq > 0, the point process x converges to \ weakly as 
m > n — > oo, m/n — >■ /3 G [l,oo) ; in which \ is a Poisson process in R + x R with 
intensity 



7T 



Ex(A x I) = — / u 2 du 

3 J A 



13 y/((l + /3- 1 / 2 ) 2 - x){x - (1 - /3-V 2 ) 2 ) 
2tt x 



-\ 4 



dx 



for any bounded Borel sets A C R+ and I e ((1 - /3~ 1/2 ) 2 + e , (1 + /3~ 1/2 ) 2 - e ). 

Remark 1.2 As in JT$, we introduce the eo here because of the technical restriction of 
the Plancherel-Rotach asymptotics of the generalized Laguerre polynomial. In Theorem 
11.11 we successfully removed the eo by using more accurate bounds. We believe the €q can 
also be removed here but it need much effort to do so. 

We can have a similar version of Corollary O Fix e > 0. Let t£° < . . . < t£° be 
the k smallest eigenvalue gaps of the form — Aj such that Aj G ((1 — /3 -1 / 2 ) 2 + e , (1 + 
0-1/2)2 _ eo )_ Moreover let 



T, 



(n) 



7T 



/V((l + /3- 1 / 2 ) 2 - x){x - (1 - /j-V 2 ) 2 ) 
27ra; 



n 1 



dx 



1/3 



(n) 



1,2..., fc. 



We have the following corollary which describes the limiting distribution of r k 
Corollary 1.2 For any < x\ < y\ < . . . < x^ < y^, then as n — )■ oo 



(«) 



(x< < n 4/3 r/ n) <y e ,l<£<k) 



e -*i - e -wS 



fc-i 



In particular, the k-th smallest normalized gap converges in distribution to with density 

P(r fc G dx) oc x 3k ~ l e~ xZ dx. 

1.3 The Unitary Universal Ensemble 

In this subsection we define the unitary universal ensemble (UUE), which is a general- 
ization of the GUE. We denote V(x) as a potential function satisfying the following two 
conditions: V(x) is real analytical on R, and has sufficient growth at x = ±oo 



V{x) 



lim 

\x\-+oo log(x 2 + 1) 



+ 00. 



In particular polynomials of the form V(x) = Yl!j=o a i x2 ^ satisfy this property. We say 
a matrix X n G C n has the universal unitary ensemble, if X n is hermitian with the joint 
probability distribution function 
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where {Aj}™ =1 are the eigenvalues of X n . Using Weyl integration, we know that the joint 
probability density function for the eigenvalues is 



p(A 1 ,...,A n )ocni^-A J | 2 e-™^ 



i<j 

This ensemble is invariant under unitary transformations (that is, p(X n ) = p(U*X n U n ) 
for any fixed unitary matrices U n ). That's the reason why we call it "universal unitary". 
Note that if we take V(x) = x 2 , then it is reduced to the Gaussian unitary ensemble case. 
But the Wishart ensemble does not belong to UUE because the potential function is only 
defined on M + 

For any potential V(x), there's a related equilibrium measure. We briefly define it 
here for integrity, but for further reference please consult [14] . Consider the minimization 
problem 



E v = inf | J J log \s - t\~ x dn{s)dn{t) + J V(t)dfi(t) 



The infimum is achieved uniquely at some measure [i = /zy, which is defined to be the 
equilibrium measure. It is proved that the equilibrium measure has a compact support 
and has form 

dfiy(x) = ty(x)dx. 

where is strictly positive in the interior of the support. 
Now we can establish the following result. 

Theorem 1.4 Suppose the equilibrium measure ^(x)dx is supported in a single interval 
[a, b] . For any fixed e > ; define 

n 

x (n) _ ^j5 (n _ 4/3|Ai+i _ Ai|jAj) ]l A . 6(o+e0i6 _ eo) . 
3=1 

Then the point process converges to x weakly as n — > oo ; in which x is a Poisson 
process in M + x K with intensity 

E X (A x J) = \ I u 2 du [ ^(xfdx 

3 J A J I 

for any bounded Borel sets A C M + and I 6 (a + eo, b — eo). 

Once again we can derive the k smallest eigenvalue gaps from Theorem 11.41 Fix 
e > 0. Let t^ < . . . < t^ be the k smallest eigenvalue gaps of the form X i+ i — Aj such 
that Aj G (a + e , b — e ). Moreover let 

r W = |^ ^(^j 1 V > * = 1. 2 ■ • ■ , k. 
We have the following corollary which describes the limiting distribution of 
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Corollary 1.3 For any < X\ < y\ < . . . < Xk < yu, then as n — >■ oo 



F(x £ < n^rf ] <ye,l<£<k) — > ( e~< - e< ) \{{yl - x 3 e ). 

In particular, the k-th smallest normalized gap converges in distribution to with density 

P(r fc e dx) oc x zk ^e~ x3 dx. 

2 Outline of the Proof 

Now we denote x^ as denned in (11.21) . As we will observe from Proposition 13. 1\ in order 
to prove converges to a Poisson point process x-, it remains to prove, for any bounded 
Borel set A C R+ and / C C, we have 

X (n \A xI)-> X {AxI) weakly. (2.1) 

However, even proving (12. ip is not easy. Following the ideas by Soshnikov in [5J, we 
consider the s-modified point process. That is, we define 

where the we only keep the eigenvalues Aj for which there is exactly one eigenvalue Al- 
lying in Xi + n~' i/4 B (recall B = {u 6 C : \u\ G A, u >z 0}). 

Here's the rationale for jf. Denote = x {n) {A x I) and N l = ^ (n) (J). Then, N 2 ^ N 
implies that at least three eigenvalues are clustered together, which is intuitively a rare 
event. Hence we should have x {A x I) — ^ n \l) asymptotically almost surely. This is 
proved in Lemma l3~3l Now in order to prove [27TI we only need to prove ^ n \l) — > x(A x J) 
weakly. This is a simplification, because rather than analyzing a two dimensional process 
X^ n \ we can turn to work with a one dimensional process instead. 

The process £ may seem complicated, but its k point correlation function pfe(Ai, • • • , A&) 
can be expressed explicitly using the inclusion-exclusion argument, see (I3.13p . Again, by 
Proposition 13. II the problem is reduced to prove that the corresponding factorial moment 
matches each other, i.e., for any integer k > 0, as n — > oo 



E 



£<»)(/)! 



(£<">(/)-*)! 



p fc (Ai, . . . , A fc )c/Ai . . . dX k -)• (J^\u\ 2 dv?J ■ Y[ l] Ai |<i. (2.2) 



Indeed, the point-wise convergence of pk to the integrand of the right hand side of (12. 2p 
is proved in Lemma 13.41 Now the only thing to do is to swap the limit and integration. 
To do this, we prove that pk can be uniformly bounded in the majority of the region (see 
Lemma 13. 5p so that we can use the Lebesgue's dominated convergence theorem. On the 
other hand, although the pk will explode in the rest of the minor region, we prove that 
since the minor region is so small, the integration of pk on that will have a negligible 
effect (see Lemma 13.61) . Combining all these, we finished our proof. 
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3 Proofs for the Ginibre Ensemble Case 

In this section we prove Theorem 11.11 and Corollary 11.11 We know from [5] that for X n 
being the Ginibre ensemble, the eigenvalues of X n j ^fn has joint probability distribution 

p(Ai, . . . , A n ) = C n exp( -n^T \Xi\ 2 J | — A ^ | 2 

^ i=l ' l<i<j<n 

where C n is a constant. Moreover we can define its /c-point correlation function as 

n\ f°° f°° 

Pkv^ii ■ ■ ■ i Afc) := -, 7TT / ••• / p(Ai, • • • , A n )c?Afc + i . . . dA n . 

(n — k)\ j_ OQ j_ OQ 

Intuitively, pk(Xi, ■ • ■ , Afc)<iAi . . . dXk is the expectation of the number of k tuples of eigen- 
values in (Ai + dXi) U . . . U (Afc + dXk). Note pk is not a probability density function as 
the total integration is not unity. The following property is quite useful in this paper. If 
we denote N as the number of eigenvalues in some interval /, then for any integer k < n 
we have 

\ n\ i r 

E (jy - fcj! = J k Pk(^i,---,^k)dXi,...,X k . (3.1) 

For a more detailed explanation, please refer to [6]. 

By using the technique of orthogonal polynomials, we can write pk in a determinantal 
form. 



p fc (Ai,...,A fc ) =n k n fe exp^-n^|A j | 2 ^det(ir„(A i A*)) 



l<i,j<k- 



Here 

n— 1 o p 

Before we begin proving the result, we first analyze the behavior of K n (x). This will 
play an essential role in later proofs. In short, for x being real, e~ nx K(x) — > 1 if x E (0, 1) 
and e~ nx K(x) — > for x > 1. For x being complex, the following lemma gives a precise 
bound on the decay rate. 

Lemma 3.1 Let z = z n = x + n~ 3 ^c. Here x E R and c E C are independent of n. 
Moreover let 

n—l £ £ 
£=0 C ' 

TTien t/ie following inequalities are valid. 

1. If\z\ < 0.02, then 

\R n (z)\ < a/^0.06". 

V 2,71 
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2. 7/0.01 < \z\ < 1, then there exists constants C and N (independent of z and n) 
such that for all n > N, 



\R n (z)\ < C^ljzle^y. 



3- If \z\ > 1, then there exists constants C and N (independent of z and n) such that 
for all n > N, 

\1-R n (z)\ < C^le-^l- 1 ^". 

Proof. Part 1. 

By applying Taylor's theorem in the complex setting on f(z) := e nz we have 

rz < _ t x„-l n e -n(z-t) 

= L (n . - Di — *■ 

Here the integral is along any curve connecting and z. In particular, we can set the curve 
to be the straight line between the two points. By Stirling's theorem, n! = v / 2vm(n/e) n e /3n 
for some (3 n G (12n + l)" 1 , (\2n)- 1 ). Thus 




a y e n : I t"- l ( - "'dt 



% - x e- ntz dt. 

V Jo 

Hence 




\Rjz)\ < xl^\z\ n e n I t n - x e- ntu dt 



in n f 

' Jo 



2tt 1 

Here we denote u = and v = If \z\ < 0.02, then e~ ntu < e n|z| < e a02n . Hence 

l*MI S ^'"^r < ^(0.02e-)"< 

Part 2. If 0.01 < \z\ < 1, then for sufficiently large n, u = U(z) > 0.005 > 0. Thus 
from part (1) of the lemma 

where g(t) = t n ~ 1 e~ ntu . By standard analysis, g(t) achieves its maximum at t* = (n — 
\)/{nu). Thus, if t* < 1, i.e., n — 1 < nu, 

\R n (z)\ < ,f^\z\ n e n g(r) = ,f^\z\ n ex V (l + (n-l)\og n '"" ' 



2tt V 2vr \ nu 

< x f^ e \z\ n < x [^ e \z\ n e^-\ z \\ 
~ V 2vr 1 1 ~ V 2vr 1 1 
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On the other hand, if t* > 1, then 



\R n {z)\ < J^\z\ n e n g(l) = J±\ z \- e ^ 




J^_\ z \n n{l-\z\) . n{\z\-u) 

2tt 1 



Let z = \z\e ld . Then exp ^n(|z| — u) j = exp ^n|z|(l — cos#) j . Since 

9 = arcsin(3r(z)/|z|) = 0(n- 3 / 4 /|z|), we obtain n|z|(l-cos0) ~ n\z\6 2 /2 ~ C^ 1 / 2 /^]). 

Hence exp fn(|;z| — w) j — )■ 1 as n — > oo. 

Thus we have complete proving the existence of these constants N and C. 

Part 3. If \z\ > 1, the quantity |1 — R n (z)\ has the following upper bound (recall 



u = ^(z)). 

n-l 



\1-Rn,(z) 



-nz 



• -E 



n— 1 < 

n r 



£=0 



n z 



Let = be the £-th term in the above summation. Then 

. ^ n-l n-2 

A? < Ai-i := CfVi 

n n n 

Thus |1 — i? n (a;)| < e~ nn A ri _ 1 ■ Yle=o ° l - Our nex ^ target is to prove Y^=o c i/V^ * s 
bounded above. If this can be done, then by Stirling's formula 

\1-Rn(z)\ < 0(l)-e- nu — JfL . ^ < 0(1) ■ (Izle 1 -") 

ni V / 

= 0(1) • (|«|e 1 - | * l ) n -e n(l * | - tt) . 

By part (b) of the lemma, exp fn(|z| — u) \ — > 1 as n — > oo. This proves the existence of 

the constants N and C. 

The only thing left is to prove ^ o q/0i is bounded above. Indeed, for £ > 
n - y/n, a < 1. For £ < n - y/n, Q < (1 - 1 / \/n) n ~^ K ~ e . Thus we get 

n-l OO y 1 \ ^ 

J> = E E ^< v^ + El 1 -^) ^ 2 v^- 

□ 

Remark 3.1 i/ie end of the proof of part 2 of Lemma A3. we actually proved a 
quite useful result. That is, e n ^^~^ — > 1 uniformly if £ is bounded away from zero and 
$*(£) = C(n~ 3 / 4 ). This will be used for several times later on. 

Remark 3.2 Since the one-point correlation function is 

Pi(Ai) = un' 1 exp(-n|Ai| 2 )if n (|Ai| 2 ), 

From Lemma / TO]) we know that pi (Ai)/n converges to 7r _1 i/|Ai| < 1 and zero if |Ai| > 1. 
27ms m agreement with the famous circular law. 
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In order to prove the convergence of the random point process % > we & Ye the 
following proposition. This is a very slight modification of the Proposition 2.1 in pQ. The 
proofs can also be found there. 

Proposition 3.1 Let x^ — Yli=i^Xi n be a sequence of point processes on IR + x C, 
and x & Poisson point process on M. + x C with intensity /i having no atoms and a -finite. 
Assume that for any bounded intervals A and I and all possible k > 1 

lim e( —. Ct^t^ttt 1 = »{A x If. (3.2) 

n->oo \(x ( - n) {A X I) - k)\J 

Then the sequence of point processes x^ converges in distribution to x- 

Moreover the following lemma is frequently used. The proofs can be found in [9]. 

Lemma 3.2 Let M be an n x n Hermitian positive definite matrix. For any u C 
{1, 2, . . . , n}, let M w (resp. M^) be the submatrix of M using rows and columns numbered 
in u (resp. {1, . . . ,n}\u). Then 

det(M) < det(M w ) det (M u ). 

Now we begin to prove Theorem ll.il that the point process x^ — J2i ^n- 3 / 4 |A i *-A i |,A i 
converges to a Poisson process. By Proposition 13.11 we just need to verify (13. 2p . To do 
this, it remains to prove that for any fixed bounded sets A C M + and JcC, 

X {n \A x /) A Poisson(A) (3.3) 

for A = 4r(/ fl \u\ 2 du)(f InD{0A) dv). 

Let A n = n~ 3 ^A = {n" 3 / 4 a : a G A}. B n = n~ 3 ^B (recall B = {u G C : \u\ G 
A,u y 0}). Consider the point process 

n-l 
i=l 

and its thinning £( n ) obtained from £( n ) by only keeping the eigenvalues A& for which 
£^(Afc + Bn) — L This technique is introduced in [5] by Soshinikov as the s -modified 
process. Intuitively, for any bounded sets A C K + and I C 1R 2 , we have 

X (n) (^x/) = #{z : Ai G /, |Ai* - Aj| G A n } 

£ (n) (/) = : Aj G J, there exists only one Xj G Aj + 

Then {%w(y4 x I) / £^(-0} would imply that there exist at least two eigenvalues 
clustering around the third eigenvalue Aj, which is a rare event. Hence we expect 
X^ n '(A x I) = £("■)(/) asymptotically almost surely. This gives us an advantage: in- 
stead of analyzing the 2-D point process x , we can turn to work with the 1-D process 
£( n ) ; which is much simpler. 

To formalize this, the following lemma says that x {n K A x J ) and £ (n) (-0 

are asymp- 
totically the same. 
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Lemma 3.3 For any bounded sets A C M + and I C R 2 , 

X (n \A x I) -l (n \l) A 0. 

Proof. Let c be such that A C (0, c) and define c n = cn -3 / 4 . 
First we prove that if 

then ^"■'(Aj + D + (0,c„)) > 2, that is, at least 3 eigenvalues are clustered together. Here 
D + (0, c n ) = {z : \z\ < c n , z >z 0} is the half disk in the complex plane. 

Indeed, if Xi* — Aj G -B n (recall the Xi* is the closest eigenvalue to Aj) and £ ( - Tl * ) (Aj+-B n ) 7^ 
1, then ^ n \\i + B n ) > 1. Thus £< n )(Aj + D+(0,c n )) > ^(A; + S n ) > 2. Conversely, 
if £ {n) (A 4 + B n ) = 1 then there exist j such that Xj - Xi G B n C D+(0,c n ). Thus 
A;. G A, + D+(0,c n ). If Ai* £ 5 n then A,* ^ A,. This implies that ^ n )(A, + D+(0, c n )) > 2. 
In summary, no matter what case happens, we always have ^ n - l (Aj + D + (0, c n )) > 2. Thus 
the statement is correct. 

From the statement we have 



rc-l 



\ X ^\A X /) -£(«>(/) | < E%")(D+(A,c„))>2 } < 5®(B). 



i=l 



where 



and 



= J2 5 (A n ,..,A Ifc ) 

Aj-p...,Aj, arc all distinct 



B = {(A, Xi, x 2 ) : A G / and (xi,x 2 ) G D + (0, c n ) 2 }. 
Since L 1 convergence implies weak convergence, we just need to prove E(S^ 3 ^(i3)) — > as 



n — > 00. We have 

E(S (3) (£)) = / dX / p 3 (A,xi,x 2 )rfxic/x 2 (3.4) 

<// </D+(A,c n ) 2 

where p 3 is the 3-point correlation function. That is, 

p 3 (A,£i,x 2 ) = n 3 7r" 3 exp(-n(|A| 2 + |xi| 2 + |x 2 | 2 ))x 

K n (\X\ 2 ) K n {Xx\) K n (Xx* 2 ) 
det I K n ( Xl X*) K n (\ Xl \ 2 ) K n (x lX *) I . (3.5) 



K n (x 2 X*) K n (x 2 x\) K n (\x 2 



\2\ 



In order to prove that (13.41) approaches to zero, we need to decompose it into three 
parts because of the different asymptotics of K n {z) in different regions. More precisely, 
we define 



h = {AG/ 
h = {AG/ 
/ 3 = {AG/ 



| A | < l- n - u ' U2 }, 

| A | > l + n-°- 02 }, 
l_ n -o.o 2< | A | <1+n -o.o 2} 



as the inner part, outer part and the middle part of /. Below we shall prove that the 
integral is small on all the three regions. Our strategy is as follows. 
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1. For the inner part, K n (z) ~ e nz . We expect that the terms in the determinant in 
(13.51) would cancel out, leaving us a small integrand. 

2. For the outer part, K n (z) ~ 0. Hence each term in the determinant is very small. 

3. For the boundary part, we can only obtain a poor upper bound for ps(Ai, A 2 , A 3 ). 
However this time the integral region of the boundary part is small, which yields a 
small result. 

Below is our rigorous proof. 
(1) Inner Part. Let's first start with I x . Replace K n (z) with e nz (l - R n (z)) in 
3.5p . we get 



p 3 (A,Xi,x 2 ) = Q n (X,x 1 ,x 2 ) + 



n 3 7r~ 3 exp(-ri(|A| 2 + |a;i| 2 + \x 2 \ 2 )) 



exp(n|A| 2 ) exp(nXxl) exp(nAa:2) 
exp(nxiA*) exp(n|xi| 2 ) expfnxia^) 
exp(nx2A*) exp^a^i) exp(n|x2| 2 ) 



(3.6) 



Here Q n (X,xi,x 2 ) is the residual, being the summation of several terms. One typical 
term is 

n 3 n~ 3 exp ^— n(|A| 2 + |xi| 2 + |x 2 | 2 ) j • exp(nAx* + nx\X 2 + nx 2 X*)R n {Xx\). (3.7) 

Taking the absolute values, we get | (!3.7p | < n 3 n^ 3 \R n (Xxl)\. Before providing a bound 
on the main part of (13. 6p . we first prove that the residual term (13.71) is small. 

First we consider the case where |A| < 0.12 is very small. In this case for n suffi- 
ciently large |Aa^| < 0.02. Thus by part (1) of Lemma 13.1 \ R n (Xx*l) uniformly decays 
exponentially, hence the impact of (13. 7p is negligible as n tends to infinity. 

Next we consider the case |A| > 0.12. Since |A| < 1 — n~ om and |xi — A| < n _3//4 , 
we obtain that for n sufficiently large, 0.01 < |xiA*| < 1 — n~ om . Thus by part (2) of 
Lemma I3.1[ the upper bound for the absolute value is 

n 3 7i- 3 \R n {Xxl)\ < 0(1) ■n 3 (|x 1 A*|e 1 - |xiA * l ) n 

= 0(1) • exp (n log |xiA*| + n — n\x\\* \ + 31ognj 

< 0(1) • exp(n log(l - n-°- 01 ) + n ' 99 + 3 logn) 

= 0(l).exp(-V 98 + o(n a98 ))^O. 

The other terms converge to zero uniformly as well. So we get that Q n (X, Xi,x 2 ) converge 
to zero uniformly — their impact is negligible. 

Next we turn to analyze the main part of (13. 6p . Expanding the determinant, we have 

p 3 (A, x u x 2 ) = vr- 3 n 3 (3 - e" n|:C2 - A|2 - e^ 1 ^ 2 - e^"^) 
+7r _3 n 3 (e™ ( ' :ClA * +X2:r i +A:!: 2~l A l 2 ~l :Cl l 2 ~l :c 2l) _ 1) 

+7r -3 n 3^ e n(xi4+x 2 A*+A4H^| 2 -ki| 2 -|z|l) — ]_) + o(l). 
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In the expression above, the first term is 0(n 3 ) x 0(ri -1 / 2 ) = 0(n 5 ^ 2 ). To analyze the 
second term, let x\ = A + n^ 3 ^ui, x 2 = A + n~ 3 ^u 2 . Then, after some calculation, the 
second term is 



7r 3 n 3 



exp^-n- 1/2 (|ui| 2 + \u 2 \ 2 -«2«I)J-1 



0(n 



5/2 N 



The same is true for the third term. Thus p%(X, Xi, x 2 ) = 0(n 5 ^ 2 ). But the integration 
domain I\ x D + (Aj, c„) 2 is of order 0(n~ 3 ). Hence we successfully proved that uniformly 



lim dX P3(A, xi, X2)dx\dx2 — 0. 



h JD+(\,c n ) 



(2) Outer Part. Next we consider the integration on I 2 . The integration domain 
is (D(n~ 3 ), but we prove that the integrand p 3 (A, Xi, x 2 ) is o(l). Indeed, expanding the 
determinant of p 3 (A,xi,x 2 ) in (13. 5p . we get a summation of several terms. One typical 
term is 

72 3 vr~ 3 exp(-n(|A| 2 + |xi| 2 + |x 2 | 2 )) • K n (\x* 1 )K n (x 1 x* 2 )K n (x 2 X*). (3.8) 

For |A| > 1 + n~ om and \x{ — A| < n~ 3 ^, we have that for n sufficiently large, l^i^l > 
1 + n~ om and \xiX*\ > 1 + n~ om for i — 1,2. Hence by part (3) of Lemma 13.11 the 
magnitude of (13.81) has upper bound 

< 0(1) -exp(3nlog(l + n-°- 01 ) - 3n a " + 3 log n) 
= exp(-3n a98 /2 + o(n - 98 )) ->■ 0. 

Similarly, the other terms in the expansion of the determinant also uniformly converge 
to zero. Hence p(X,Xi,x 2 ) = o(l) uniformly. Thus we conclude that the integral on I 2 
also converge to zero. 

(3) Middle Part. Finally, we prove the integral on I 3 is also of o(l). Again we 
expand the determinant and prove that each term is small. Since for any zeC, 

n—1 p \ if oo p I I p 

iA-„wi<£^<£^ = e»'«', p-9) 

£=0 £=0 

Expanding the determint in (13.51) again and using the estimate (13.91) . we have that the 
magnitude of one typical term in the expansion has upper bound 

n 3 7r~ 3 exp(-n(|A| 2 + |xi| 2 + \x 2 \ 2 ))\K n {Xx\)K n (xix* 2 )K n (x 2 X*)\ 
< n^expj-^ 1 ~ |Xl|)2 - " (|A| ~ |X2|)2 - ~ |X2|)2 U Q(n 3 ). 

But the integration domain is of order 0(n~ 3 ) x O(n~ ' 02 ) = 0(n~ 3m ). Hence the total 
integration over J3 decay to zero at the speed of O{n~ om ). 

From all the above, we proved that the integral of I = I\ U I 2 U I3 converges to zero. 
Hence E(E^(B)) -> 0. □ 
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Note that previously our ultimate goal is to prove y^iA x I) — >■ Poisson(A) for 
A = l 701 ^ ' 1 )! J A r 3 dr. Now Lemma 13731 tells us x^ n \A x I) and £( n ) has asymptotically the 
same distribution, it remains to prove ^ n \l) — > Poisson(A). Or, we just need to prove 
that 



E 



£(»)(/)! 



L(eW(/)-*)! 



— ^ |w| 2 cfii 



/nD(o,i) 



du ,Vfc > 1. 



(3.10) 



This is because if (I3.10P holds, then all the moments of converges to the moment 

of Poisson(A). This implies ^ n \l) — > x(A x /) weakly, or x^ n \A x /) — y x(A x /) weakly. 

Denote Pfc(Ai, . . . , A&) as the /c-point correlation function of ^ n '. Recall the point 
correlation function has the property (13.1 1) , the left hand side of (I3.10p can be expressed 
as _ 

£<»)(/)! 



E 

■.(£(»>(!)-*): 

Thus now our ultimate goal is to prove 



Pfc(Ai, . . . , AjtjdAi . . . d\k- 



(3.11) 



Pfe(Ai, . . • , Xk)d\\ . . . d\k 



1 



7T 



\u\ du 



B 



dv) ,Vfc > 1. 



/nD(o,i) 

To establish (I3.12p . it remains to prove the following three statements. 
1. If all the Afc's are distinct and have magnitude not equal to 1, then p(Ai, 

(^Lkl 2 ^) fc -n- =1 %i<i- 



(3.12) 



, A fc ) -> 



2. Define 
Then p(Ai, 



{(A 1; . . . , A fc ) G J fe : (Aj + B n ) n (A, + B n ) = 0, 1 < i, j < fc}. 
. , Afc) is uniformly bounded in f2 (Lemma 13.51) . 



3. Define f2 = 1 then the contribution of the integral of Pfc(Ai, . . . , Aj.) in is 
negligible since the volume of f2 is sufficiently small (Lemma I3.6p . 

As a little discussion, the set Q represents the majority of the integration region. 
Combining the point-wise convergence result in statement 1 and the boundedness result 
in statement 2, we can perform the dominated convergence theorem in the main region 
Q. Furthermore statement 3 ensures that the contribution of the minor region Q does 
not play an important role, hence we can obtain the result (I3.12p . 

We prove the three statements one by one in the following. 

Lemma 3.4 (Point-wise convergence) Let Ai, . . . , A& be distinct complex numbers and 
| Aj | 7^ 1 for i = 1, . . . , k. Then as n — >■ oo 



Pfc(Ai, . . . , Afc) > — r 



7T 



B 



k k 

\a\-da ) ]Jl 
i=i 



|Ai|<l- 
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Proof. To begin with we consider the first case where |A;| < 1 for all i. Since all the Aj's 
are distinct, for n large enough, A« + B n are disjoint sets for different z's. Moreover there 
exists some e > such that for n sufficiently large, Aj + B n C D(0, 1 — e ). 

By an inclusion-exclusion argument, for sufficiently large n we can represent pk in 
terms of pk as follows. For a general result please see [5]. 

Pk(Xi, ■ ■ ■ , A fc ) = ^ — / dxi... dx k 

m =0 J\i+B n J\ k +B n 

I P2k+m{Xi,xi,...,Xk,x k ,yi,...,y m )dy 1 ...dy m . (3.13) 

J ((Ai+B n )U...U(A fc +B n )) m 

Since p2k+m = for 2k + m > n, (13 . 13[) is a finite sum thus there is no convergence 
issue. First of all we consider the m = case. Now the expression for p 2 fc is 



P2fc(Ai, X\, . . . , Afc, Xfc) 



n 2fc 7r 2/c exp( -7i^(| A 



■" ' 'l<<J<* V K n(Xi\*) K n (XiX*) 



,2k -2k 



" 7 i<S<* V e-"(l^l 2+ l^l 2 )/ 2 K n (x,A*) e-^l^l 2 +l^l 2 )/ 2 K n (x^*) 



(3.14) 



Here the determinant is (2fc) x (2k) with sub 2x2 blocks described above. We prove 
next that only the diagonal 2x2 blocks can have a non-negligible contribution. Indeed 
consider the (i, j)-th 2x2 block. If i ^ j, then it has the estimation 



e ^M 2 M*j\ a )ftK n (xi\f) e-^^+^^KnixiX*) 



0(1) (1) 
0(1) 0(1) 



This is because, for example consider the top-left entry, 

e -n(\^ + \^ 2 )/2 Kn{Xi X*) ~ exp(-n(|A l | 2 + |A/ - 2A*A i )/2) = o(l) 

as Xi Xj. Moreover the o(l) notation decays to zero exponentially fast. The same is 
true for other terms. Hence the determinant of the small block above tends to zero. 
However, if % = j, then 



e- n ^ 2 K n (\X. t \ 2 ) e- n ^ 2+ M 2 y 2 K n (XiX*) 

e -n(\x^+\\^)/2 Kn ( x .X*) e-^KnQXil 2 ) 



1-e 



-n\\i-Xi\ 2 



This is of order 0(n~ x l 2 ) compared to the exponentially decay in the i ^ j case. As a 
consequence, in the expansion of the determinant in (I3.14p over all permutations of S2k, 
only the terms consists of the entries in the diagonal 2x2 blocks can make a non-trivial 
contribution. Indeed, their contribution is exactly 



JJ / p 2 {Xi,Xi)dxi ~ Yl 

1=1 J^i+B n i= l 



2 -2 
U 7T 



(1 — exp(— n\Xi — Xi\ 2 ))dx t 



\+B n 
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(3.15) 



Let Xi = Aj + n 3//4 Mj, where U{ G B, then 

(j335D = (n 1/2 TT- 2 J (l-exp(-n- 1 / 2 | M | 2 ))^ ~ ( 7r_2 / l M l 2rfM ) • 

Having analyzed the m = case, next we prove that the contribution of the terms 
corresponding to m > 1 is negligible. 
By lemma [3T21 we have, for m > 1, 

m 

P2fc+m(Ai, a?i, . . . , A fe , x fc , yi, ... , y TO ) < p 2 fc(Ai, £i, . . . , X k , x k ) pi(yi). 

8=1 

Thus the contribution for m > 1 is bounded by 

/ dxi... I dx k p 2k (Xi, a?i, - - - , A fc , x k ) ) ^ — ( / pi(y)d?/] . (3.16) 
J\ 1 +B n J\ k +B n ' m >\ m [ \Jb n J 



Here 6 n = (Ai+£> n )l_l. . .U(Afc+-B n ) has size C(n 3 / 2 ). However since |i^ n (z)| < exp(n|z|), 
we have 

\ Pl (y)\ = mr' 1 exp(-n\y\ 2 )\K n (\y\ 2 )\ < mT 1 = 0(n). 

Thus the second factor of (13.161) will tend to zero. Since the first factor is just the 
m = case, which converges as proved. Thus the whole expression converge to zero. 
Combining our result for the m = and the m > 1 case, we successfully proved that 
p fc (Ai, . . . , Afc) ->• J B \u\ 2 du) k when |A;| < 1. 

As the second step, we prove that for distinct Aj's, if there exists some io such that 
|A.; | > 1, then pfc(Ai, • • • , Afc) — > 0. Indeed, we can assume that Ai, A2, . . . , A p > 1 + eo 
and A p+ i, . . . , Afc < 1 — e for some e > and p < k. Then from the first part, we have 
shown in (I3.16P that 

Pfc(Ai,...,Afc) < ( / dx x ... I dx k p 2k (Xi,xx,...,X k ,x k ) )x 



m>0 

Again by Lemma 13.21 we have 

k 

P2k(X 1 ,x 1 , Afc,Xfc) < J^p 2 (Ai,Xi). 

i=i 



£^(7 My)dy) m . (3.17) 



Thus 

Pk 



(Ai,...,Afc)<TTf / p 2 (X i ,x)dx)y2—( [ pi(y)dy) . (3.18) 



We have already shown that for i > p, 



p 2 (Xi, x)dx — > 7r / \u\ du 

Xi+B n J B 
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which is finite. But for i < p, for n sufficiently large, for all x G Aj + B n , for n sufficiently 
large we have \x\ 2 > l + e /2. Then by Lemma [3. II we have 

|^n(M 2 )| < 0(1) ■e" |a;|2 (|a;| 2 e- (|a;|2 - 1) ) n < 0(1) ■ e n ^ ((1 + e /2) e - eo/2 ) n 

Hence for i 6 Aj + B„, 



|A,| 2 


— n 


\x\ 


|A,| 2 


— n 


\x\ 



< n 2 ir~ 2 exp(-n|A<r - ra|a;p)lif r »(|Ai| je )iif w (|arp) 

< 0(1) • n 2 ■ m + e /2)e- eo/2 J ->■ 0. 

Thus nLi J\ +b P2(K,x)dx — > and the second factor in (I3.18P is bounded as having 
been previously shown. Hence we have Pfc(Ai, . . . , Afe) — > 0. The proof is complete. □ 

Lemma 3.5 (Uniform boundedness) There is a constant C depending only on the set 
A such that, for any n > 1 and (Ai, . . . , A&) G fi, 

p fc (Ai, . . . , A fc ) < C. 

Proof. In (I3.18P we have shown that 



p 2 (Ai,x)da; J V — ( / 



Pk(M, • • • > A/t) < 

As is proved in Lemma [3.41 the second factor does not depend on A and is convergent. 
We just need to prove that the first product is uniformly bounded. We note 

n , 2-2 A J e-^X(|A| 2 ) e-^l 2+ l A l 2 )/ 2 ^(Ao;*) \ 

We consider the following three regions separately and shall prove that p2{X,x) is 
uniformly bounded on all of them: (1)|A| < 0.12, (2) |A| > 100 and (3)0.12 < |A| < 100. 

Part 1. If |A| < 0.12, then as n is sufficiently large, |A| 2 , \x\ 2 and |Ax*| are less than 
0.02. Writing K n (z) as e nz {l — R n (z)), by part (1) of Lemma l3TTj we get 

( 1 e -n(M 2 +|A| 2 -2Ax*)/2 \ 

p 2 (A,x) = n vr det I _ nQx \2 +W 2_ 2xX . )/2 l J + o(l) 



(l - exp(-n\x - A| 2 ))+o(l) 
0(n 2 ) ■ Otn- 1 ' 2 ) + o(l) = £>(n 3/2 ). 



But the integration domain is 0(n~ 3 / 2 ). Hence the integration of p2(X,x) in A + B n is 
uniformly bounded. 

Part 2. If |A| > 100, then for sufficiently large n, we can ensure \x\ 2 , |A| 2 , \x\*\ > 
100. By part (c) of the Lemma I3TTT this implies \K n (\)*\ < C ■ (lOOe -99 )™ for some 
constants C. Hence 



n 



2 n- 2 e- n ^ 2 -^ 2 \K n (\x*)\ < Cn 2 n- 2 ■ (lOOe 



99\n 
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which is uniformly bounded above. The same is true for other terms in the determinant. 
Hence the integration of p 2 (X,x) is bounded, just as in Case 1. 

Part 3. Finally, we prove the most tricky part: 0.12 < |A| < 100. For p 2 (X,x), we 

have 



p 2 (X,x) 



n 2 n 2 det 



n 



e- n l A l 2 fs: n (|A| 2 ) e~ nXx * K n (Xx*) 
e~ nxX * K n \xX*) e~ n \ x \ 2 K n \\x\ 2 ) 



(3.20) 



First we analyze the second term in (I3.20p . Using the inequality |-fT n (z)| < e n ' 2 ' we have 

n 2 7T- 2 \e- 2n ^ xX ^ -exp(-n\X\ 2 -n\x\ 2 )\ ■ \K n (Xx*)\ 2 
< n 2 7T- 2 \e- 2n ^ ixX * ) - exp(-n|A| 2 - n|x| 2 )| • exp(2n|Ax*|) 

(2n{\x\*\ - M(x\*)))-exp(-n(\\\ - |x|) 2 



^ 2-2 



expl 



, 2-2 



exp(2n(|xA*| - 9?(xA*)) 



i 2 -2 



exp^— ?7.(| A| — \x\)' 



(3.21; 



Here (I3.2ip refers to the line next to it. For the first term of (I3.2ip . we have \x\*\ — 
9£(xA*) = 0(n~ 3/2 ) (see Remark EUj) . Thus the first term is of order 0{n 2 ) ■ 0{n~ x l 2 ) = 
0(n 3 / 2 ). 

Moreover, for the second term of (I3.2ip 



2 -2 
U 7T 



fl _ e -n(|A|-W)^ = C ( n 2) . fa _ e 



_ p -0(n"V2) 



C»(n 2 ) • 0(n 



0(n 



3/2 \ 



Hence (I3.2ip = (9(n 3//2 ) uniformly. Since the integration region is of order (9(n _3//2 ), we 
obtain that the integral of the second term in f )3.20p is 0(n~ 3 ! 2 ) • 0(n 3 ^ 2 ) = 0(1) which 
is uniformly bounded. 

Now it remains to prove that the first term in f)3.20p is 0(n 3 ^ 2 ) hence is uniformly 
bounded after integration on the region of order 0(n~ 3 ^ 2 ). 

Define f(z) := e~ nz X^=o n ^ ' z& ' I • Then the determinant is nothing but 



/(|A| 2 ) /(Ax*) 
f{xX*) f(\x\ 2 ) 



fW) f(Xx*)-f(\X\ 2 ) 
f(xX*) - f(\X\ 2 ) f(\x\ 2 ) - /(Ax*) - f(xX*) + /(|A| 2 ) 



(3.22) 



the equality holds because the determinant is invariant after subtracting the first row 
from the second one and then subtracting the first column from the second one. In order 
to prove p2(A,x) = 0(n 3 ^ 2 ), we just need to prove that the determinant in (I3.22p is of 
order 0{rT x l 2 \ A simple calculation with Stirling's formula gives 



—ne 



,n n z n 1 




lIL e -^ e n(l-z) n-l 

2vr 



for some (3 G 



i i 



12n' 12rt+l 



). we have 



/(xA*)-/(|Ap 



xX* 



/'(Ode 



(3.23) 



I A! 2 
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where the integration is along the straight line connecting the two end points. For £ lying 
between |A| and xA*, 3(f) = 0(n^l 4 ). 

|f (01 < 0(1) • ^(e 1 ^^) V^H5I) < 0(1) ■ yfr- e «(*«)-lfl). 

By Remark 13. 1[ e 71 *- 3 *^ - '^ = 0(1) if f is uniformly bounded away from zero and 
3(f) = 0(7?r 3 / 4 ). This can be achieved when |A| > 0.12. Hence we have /'(£) = 0(y/n) 
uniformly. However, in (13. 23j) the length of the integration region is |A||x — A| = 0(n~ 3 ^). 
Thus we have f(xX*) - f(\X\ 2 ) = C^n" 1 / 4 ). Similarly /(Ax*) - /(|A| 2 ) = 0(n- l ' A ). As 
a conclusion, the second term of the expansion of the determinant in f)3.22p (/(Ax*) — 
/(|A| 2 ))(/(xA*) - /(|A| 2 )), is of order 0(rT x l 2 ). 

In order to prove that the first term in the expansion of (13.221) is also of order 
0{n- l l 2 ), we just need to prove /(|x| 2 ) - /(Ax*) - /(xA*) + /(|A| 2 ) = 0(rT x l 2 ) because 
/(|A| 2 ) is bounded above. Now 

/(|x| 2 )-/(Ax*)-/(xA*) + /(|A| 2 ) 

= y V/'C&O - A*/(eA*))de = f £ (f'itv) + tijffa^dndt (3.24) 



and 



,n n z n 1 r,, , N n 2 



f(z) + zf"(z) = e -nz'_^J__ . n 2( z _ X N ^^^n-l (2 _ ^(1-*) (3 35) 

n\ e^V27rn 



2h 



where ft, e (i^, i^+t)- 

Since the integration region of (13.241) is 0(n~ 3 ^ 2 ), we just need to prove that f'(z) + 
zf"(z) = 0(n) for z such that 0.01 < \z\ < 100 and 3(z) = C(n- 3 / 4 ). 

Thus by (I3.25P the only thing left to do is to prove (z — l)z n ~ l e n ^-~ z ' = 0(n~ 1 / 2 ), or 

(Z - l)2 z 2n e 2„(l- 2 ) = 0( n -l) (3 26) 

Let \z\ — 1 + u, and z — (1 + w)e . Then 

|(2- l)^ n e n(1 - 2) | 2 = 0(1) • |z- l| 2 ^(l + M)e-" N 
where we used Remark 13.11 again. Now since 

= C(n~ 3 / 4 ) and \u\ < 100, 
\z-l\ 2 = u 2 + 2(l + u)(l- cos9) <u 2 + (l + u)9 2 <u 2 + 1019 2 , 
thus \z — 1| 2 = u 2 + 0(n~ 3 / 2 ) and the constant in O(-) is uniform in z. Hence 

\(z- l)z n e n{1 - z) \ 2 < 0(l)-u 2 ((l + u)e- u y n +0(n- 3/2 ). 



Let g(u) = u(jl + u)e 11 ^ . Then it is easy to obtain that g(u) achieves its maximum 

when u = (1 + \/l + 4n)/(2n) = (9(n -1 / 2 ). The maximum value of g(w) = (9(n -1 / 2 ). We 
obtain 

\{z - l)z n e n{1 - z) \ 2 < 0(1) ■ g(uf + C(n- 3 / 2 ) < C^ 1 ). 
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This is exactly (13. 26 p . Thus we successfully proved that f"(z) = 0(n), or, the deter- 
minant in (I3.22p is of order 0(n~ l l 2 ). We conclude that the integration of p2(X,x) is 
uniformly bounded. 

From all the three parts above, Pk(X±, . . . , Xk) < C for some constant C when 
(Ai, . . . , Afc) 6 0. □ 

Lemma 3.6 (Negligible set) Let Q be the complement of Q in I k . Then asn-^oo 



Pfc(Ai, . . . , AfcJdAi . . . dAfc 



0. 



Proof. The strategy for this proof is to show that \Q\ decays sufficiently fast — much 
faster than the growth of Pfc(Ai, . . . , Xk). Indeed, if (Ai, . . . , A&) G Q, then this means there 
exists some clusters of at least three eigenvalues, which is a rare event. To formalize this 
idea, we need to introduce a rigorous definition of the cluster. 

We define a equivalence relation in {Ai, . . . , A^}. For any Aj, Xj, if |Aj — Xj\ G B n then 
we say A.; ~ Aj. If there are a series of points A, 1; . . . , X it such that Aj ~ A^, X ia ~ X is+1 for 
s — 1, 2, . . . , t — 1 and X it ~ Xj, then we also say Aj ~ Aj. In the end, we define Aj ~ Aj. 

Suppose there are p equivalent classes. From each of them we draw representatives 
At 1( . . . , Aj p , which is the largest number in that class with respect to the ordering (•-<•)• 

Geometrically, each equivalent class represents a cluster of eigenvalues. For a more 
intuitive graph, please see Figure [TJ In this figure from the 25 eigenvalues on the complex 
plane, we have 4 equivalent classes (clusters). In each of them we have a representative 
eigenvalue which is at the top of the cluster. 



1 * N 



1 X | Equivalent Classes 
\ / (Clusters) 



® 



Representative 
Eigenvalues 



Figure 1: Equivalence Classes and Representative Eigenvalues. 
Then we have the obvious bound 



p k (X h . . . , X k ) < I dxi 



< n u 



d>XpPk + p(Xi, . . . , Afc, Xi, . . . , Xp) 
P2(X i] ,x j )dx j . 



j = l " ^'j 
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where we used Lemma 13.21 in the second inequality. As is proved in the Lemma 13.51 the 
second product is uniformly bounded in Aj. and n. Thus we only need to deal with the 
first product. 

But pi(Ai) = 0{n). Thus p(Ai, . . . , A/.) = 0{n k ~ p ). But the region of the integral is 
of 0(n~ 3 ( k ~ p ^ 2 ), thus pk is of order 0{n~^ h ~ p ^). Since k > p, otherwise (Ai, . . . , A&) will 
belong to fl, we conclude that 

p fc (Ai, . . . , \k)d\i ...d\ k — > 0. 

□ 

Proof of Theorem 11.11 With Lemma [331 E31 GLS1 the proof of (I3.10p is quite 
straightforward. 

Pfe(Ai, . . . , \k)d\i . . . Afe 

loPfc(Ai, . . . , \k)d\i . . . A fc + / %p fc (Ai, . . . , \ k )d\i ...d\ k 
jk jjk 



The first term is uniformly bounded by Lemma 13.51 thus we can use the dominated con- 
vergence theorem. The second term converge to zero as stated by Lemma 13.61 The proof 
is complete. □ 

Proof of Corollary 11.11 We first denote tp < . . . < tp to be the k smallest 
elements in the set {|Aj — Aj* | : 1 < i < n}. Correspondingly we define = (7r/4) 1//4 t^ 
for 1 < £ < k. We first analyze the behavior of instead. 

The event {xi < rc'^T { f ) < ye, 1 < £ < k} is equivalent to 
x H ((4/7r) i/4 (xfc)2/fc))C) > 1; 

x (B) ((4A) 1/4 (^.w)»c) = M<*<*-i, 

X (B) ((4/T) 1/4 (y*-i,^),C) = 0,!<£<k. 
Here i/q = 0. The independence property of the Poisson process yields that 

F{x e < n 3/4 r^ n) < y t , 1 < £ < k) 



U _ e -ivt-<)\ JJ (M J _ 4 )e -(^-4) f[ e -i4~vU) = L-< - e -vt\ J[(y 

/ 0-1 V / 



; 



Let %i — > yi~, we can get that the density function for n 3 / 4 (r{ . . . , ), in the limit, 
is proportional to 



u\...u\e <. 



Thus, by integration w.r.t. u\, . . . , Uk-i we can get the distribution function for n 3 / 4 -^. 
P(^fc <x) = Cfc / u\e~ u *duk I u\. . . Ufr^dui . . . duk-i 

J0 J0<u 1 <...<u k 

= c fc I uf~ l e~ u *du k j vf...vl_ l dvi...dv k - 1 . 

JO Jo<v 1 <...<v k _ 1 <l 

23 



We see that P(r| n ' ) G dx) = x 4k 1 e xi dx. To prove that also has the same 
distribution, we just need to prove 



P(^ n) 7^ t^ n) for some K fc) as n -4 oo. 



Recall that if is the £-th smallest entry in the set {|Aj — Xj\ : % ^ j}. Now we 
denote \X Pe — X qe \ = t'f 1 . Then our first observation is that if there does not exist three 
eigenvalues A» 13 Aj 2 , Aj 3 such that |A,, 2 — AjJ, | A» 3 — A» x j < 2tjj. , then the 2k points {X Pe , X qt : 
£ = 1, . . . , k} must be distinct. This further implies that = Vf^ for £ — 1, . . . , k. In 
short, we have 

{t^ 7^ for some £ < k} =>- {3A ii; A i2 , X i:i such that \X i2 — AjJ, |Aj 3 — AjJ < 2t[ n ^} 
If we denote 

^i 1 ^%2'^ i 'j, are distinct 

and 

B = {{X,xx,x 2 ) : |A| < 2, |xi - A|, \x 2 - A| < Mn" 3/4 } 
where M is a large constant, then we must have 

P(tJ n ) ^ ^ for some £ < k) 

< P(H {3) (£) ^ 0) + P(?[ n) > Mn" 3/4 /2) + P(there exists some |A,| > 2). 

(3.27) 

< E(2 (3) (£)) + P(4 n) > Mn" 3/4 /2) + P(there exists some |A*| > 2). 

From the proof in Lemma 13.31 we know E(H^ 3 ^(i3)) — > as n — > oo. From the circular 
law we also know that P(there exists some |Aj| > 2) — > 0. Thus taking n — )■ oo in (I3.27P 
yields 

limsupP(^ n) ^ tf ] for some £ < k) < ¥(n 3/A r k > (tt/4) 1/4 M/2) 
Finally taking M — > oo proves the result. 



4 Proof for the Wishart Ensemble Case 

Next we consider the complex Wishart ensemble. For A mn ~ W2(m, n) being the Wishart 
ensemble, the joint distribution function of the eigenvalues of A mn is (see [10], [TTj ) 

p(Ai, . . . , A n ) oc H \Xi - Xj\ 2 J[ X?~ n exp I -m J>J. 

i<j i=l ^ i=l ' 

The /c-point correlation function is given by 

Pk(Xi, . . . , X k ) = ( K n (Xi, Xj) ) 

V / i,j=l 
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where 



K n {x, y) 



tff n -i(mx)ilf n (my) - ip n (mx)ip n _i {my) 



y/mn- (4.1) 

x-y 

m 3/2 n 1/2 ^tp n (my) J il)' n _ x (mzt)dt - ^^(my) J ip' n (mz t )dt^ (4.2) 



Here z t = tx + {1 — t)y and 



and Z^ m n ^ is the associated (generalized) Laguerre polynomial, i.e., 



oo 

e" 



- x x m - n L^- n \x)L^- n \x)dx = (P + m n ) l 6 



pi pr 



Now we denote ft := m/n, we consider the case where m,n — > oo but their ratio 

0e[l,oo). 

Note that in the Wishart ensemble case, thing are much simpler than the Ginibre 
case — all the eigenvalues are real. We can use exactly the same scheme as in section 2 
or in PQ. The only difference is the kernel. In the following, we will argue Lemma 13.31 — 
13.61 still holds. For notational simplicity, we denote a = (1 — /3" 1 / 2 ) 2 and b = (1 + /3 -1 / 2 ) 2 . 

Again, for any e > we define 

n 

£(») :=Y^ S Ai£(a+e ,b-e ) and e (n) = ^26 Xi t Xie(a+eo>b _ eo) 
j=l * 

where ^ n - ) is the thinning process by only keeping Aj's such that £(Aj + n~ A ^A) = 1. To 
show x {A x I) = £ (J) asymptotically almost surely, the same argument in Lemma 13.31 
still works. The only difference is to prove that for A G (a + e , b — e ) and |A — x\, |A — y\ = 

C(n- 4 / 3 ), 

/ K n {\,\) K n (X,x) K n (X,y) \ 
det K n (x,X) K n (x,x) K n {x,y) = 0(n 7 ' 3 ). (4.3) 
\ K n (y,X) K n {y,x) K n (y,y) ) 

Indeed, after subtracting the first column from the second and the third column, we 
obtain that the entries in the second and the third columns are (9(n~ 4//3 ) ■ \d x K n \. If 
we can prove \K n \ = 0(n) and Ic^i^l = 0(n 2 ), then each term in the expansion of the 
determinant of (14. 3 p is of order 0(n 7 ^ 3 ). The upper bound for K n and dK n is proved in 
Lemma 14.21 

The point-wise convergence the correlation function, which is a similar version of 
Lemma I3~4"l is the most difficult part. In Lemma I4TT1 we see that the diagonal 2x2 blocks 
can have nontrivial contributions, which is of order (9(n 4//3 ). Moreover we calculate this 
contribution explicitly. To prove that the contributions of the off-diagonal 2x2 blocks 
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are negeligible, it remains to prove \K n \ = 0(n), \dK n \ = 0(n 2 ), \d 2 K n \ = 0(n 3 ) when 
| a; — y\ = 0(n~ 4 ^ 3 ) and \K n \ = 0(1) when \x — y\ is bounded away from zero. This is 
because we can subtract the (2j — l)-th column from the 2j-th column and then subtract 
(2j — l)-th row from the (2j)-th row, for all j = 1, . . . , k. Then all the entry in the 
off-diagonal 2x2 blocks are of order 0{1) while the diagonal 2x2 blocks are of order 



\K n 



\d y K n \ ■ 0(n 



-4/31 



\d x K n \ ■ 0(n~^) \dl v K n \ ■ 0{n-W) 



0(n) 
0(n 2 ' 3 ) 



0(n 2 ' 3 ) 
0(n^ 3 ) 



(4.4) 



Hence in the expansion of the (2k) x (2k) determinant, if one term involves off-diagonal 
terms, then it is of order (9(n 4fc / 3-1 / 3 ), which is negligible. We note that the upper bounds 
of \d^K n \ are shown in Lemma [4.21 for j = 0, 1, 2. 

The uniformly boundedness result (Lemma 13. 5j) and the negligible result (Lemma 
13 .6 j) can also be obtained using exactly the same argument. The only difference is to 
prove that for x, y G (a + eo, b — cq) and \x — y\ < (9(n~ 4,/3 ), 



det 



K n (x,x) K n (x,y) 
K n (y,x) K n (y,y) 



0(n 



-4/3 n 



(4.5) 



Indeed, we use the same trick again. We subtract the first row from the second row, and 
then subtract the first column from the second one. Then (14. 5|) exactly becomes (I4.4|) . 
Hence the result holds. 



In conclusion, Theorem II .31 holds if we can prove Lemma [4. II and Lemma [4.21 below. 
Now we conquer them one by one. 

Lemma 4.1 If \x — y\ = 0(n~ A / 3 ) and x G (a + eo, b — e ), then 



det 



K n (x,x) K n (x,y) 
K n (y,x) K n (y,y) 



— >■ 



■242 
-7T n u 



(3 v /((l + /3-l/2)2_ x )( x _ (1 _ 0-1/2)2) 

2n x 



-0(n) (4.6) 



where the O notation is uniform in x, y and u := y — x . 



Proof. The proof is quite lengthy and just involves messy calculations. By the Plancherel- 
Rotach asymptotic as described in [2], we have 



£(ro-n) ( 



rax) 



-/V™ sin O 



X 



mx/2 e -n(f3-l)/2 % -n{^l) /2 x 



sin 



I 71 T) f) I 7T" 

(n + ~)0 O - - sin 20 o + -f sin29 1 - (nf3 + -)0 X + - +0( n -' 



where the angles 9 ,9i are defined by 



cos 0Q - 



P 



fix 



0o G (0,tt) 



2jfii 1 

sin0! := sin0 o //3, 0^(0,^/2). 



(4.7) 



(4.8) 
(4.9) 
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By defining (3 n ,x via 



n((3 — 1) = (n — l)(/3 — 1), nx = (n — l)x 



(4.10) 



we can obtain a similar formula for L^^ 1 ' (mx), stated below. 



- (m— n) 



1 



mx) 



- (m—n) . 
J n-\ 1 

. e nx/2 e -nC8-l)/2gn/2--n(^-l)/2 x 



sin 



n 



x 1 / 4 A/7rnsin #o 

It? t? n I 7T " 

N - - sin 2^ + — sin 2^ - (n/3 - -)Q X + - +C(^- 1 ; 



_ _ , (4.11) 

2' 2 2 v ' 2' 4. 1 

For notational simplicity, we define the angle inside the sine function of (14. 7p . (14. lip 
to be M + (0 o ,i) and M_(0 o ,i), respectively. Finally we can plug (JUTJ) and f fl~TT]) in f l4~T|) 

to get 



K n {x,y) 



1 



a; — y) (a;?/) 1 / 4 7r -^/sin #0 sin <p Q 
sinM_(^ o ,i)sinM + (0o,i) - sinM+(# ,i) sinM_(0 o ,i) + 0(n~ l ) 



(4.12) 



Here the angles 0o,0i are the counterparts of (14. 8p . (14. 9 p with respect to y. Moreover, 
using the same technique in [3] by using higher order expansions, we can further show 
that the error term in (I4.12p is 0{x — y) (this is intuitive because the error term is a 
symmetrical difference with respect to x and y). 
By the equality 

sin(a — b) sin(c + d) — sin(c — d) sin(a + b) = sin(a + c) sin(<i — b) + sin(a — c) sin(6 + d) 
we have 



K n {x,y) 



(x — y) (xy) 1 / 4 7Ta/ sin 6* sin O 



C(x — y) — sin 



?o -t- <PoJ 



7 i -r </>i. 



sm 



^(^o-0o)--(sin2^o-sin20 o ) + — (sin2^ 1 -sin20 1 )-n/3(^ 1 -0 1 ) \. (4.13) 



Let y = x + m where w ~ 0(n 4//3 ). Then we can expand 0o> 0i near 9o, We list 
the result below. 

\//? COS fli 2 

00-^0 = -z — —z-u + Oiu), 
2x sm 6*o 

COS 6*0 sn, 9\ 

We then plug these approximations in (I4.13P to get 

K n {x,y) 



u (xy) 1 / 4 7r y/ sin 6 sin <p 
1 



(xy) 1 / 4 7r ^sin 9 sin O 
1 



\/x7r sin 6q 
1 



C(m) + TCy//3xg(x) sin irnug(x) + 0(nu 2 ) 

Oil) + ixJ~(3xg(x) irng(x) - ^n 3 n 3 u 2 g 3 (x) 
L 6 

1 



0(1) + 7r\//3xg(x) irngix) - ^7r 3 n 3 u 2 g 3 (x) 

6 



n#(:r) - -7r 2 n 3 u 2 g 3 (x) + 0(1). 
6 
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where g(x) is the density for Marceko-Pastur Law 



(3 v /((l+/3-V2)2_ x ) (a _ (1 _ 0-1/2)2) 



9K X '■= 7T ' 

27T X 

Finally, by using the approximation, we can obtain the determinant 

K n (x,x)K n (y,y) - K n (x,y) 2 = ^n 2 g 4: (x)n 4 u 2 + 0(n). 
This completes the proof. □ 

Lemma 4.2 Uniformly for x,y G (a + eo, & — eo), K n (x, y) = 0(n). Moreover, the first 
and second partial derivatives of K n (x,y) are of order 0(n 2 ) and 0(n 3 ), respectively. 
Finally, if in addition \x — y\ > 5 for some constant 5, then K n (x,y) = 0(1). 

Proof. By the Plancherel-Rotach asymptotics for Laguerre polynomials (14. 8 p and the 
definition for ipi(x), we obtain that ip n -k(nx) = 0(n~ l l 2 ) for fixed integer k. By the 
formula 



d . , , rn - n + 21 , . , J 1(1 + m-n) . . . 
-^ipe( x ) = ~2^ l \ x ) + 2x x ^-H x )' 

we obtain that vp'(n — k)(nx) = 0(n~ x l 2 \ Taking derivatives again yields ijj"(n — k)(nx) = 
0(n~ l l 2 ). If we plug this in (14. 2p . we observe that K n (x,y) = 0(n). Taking partial 
derivatives of K n (x,y) gives that the first order derivative of K n (x,y) is of order 0(n 2 ) 
and the second order is of 0(n 3 ). 

On the other hand, if \x — y\ > 5 is uniformly bounded from zero, then from (14. ip 
we conclude K n (x,y) = 0(1). □ 



5 Proof for the Universal Unitary Ensemble Case 

For the UUE case, we proceed by using exactly the same argument as in the previous 
section. The fc-point correlation function has a similar formula 

Pk = (K n (X h Xj)) 

V / i,j=l 

where the kernel K n (x,y) can also be defined via orthogonal polynomials. 

n-l 

K n (x,y) = e-^ v W +v ^J2Pi( x )Pi(y)- 

3=0 

Here pj (x) is the j-th orthonormal polynomial with respect to the weight e~ nV ^ . 

Throughout the whole section we assume that the equilibrium measure ^(x)dx is 
supported on a single interval [a, b}. 

To prove Theorem 11.41 we need to run the argument again. As is analyzed in the 
previous section, we just need to prove the following two Lemmas. 
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Lemma 5.1 For any eo > ; if \x — y\ = 0(n 4 / 3 ) and x G (a + eo, & — eo), then 

/ Jjf n (z,z) ^,!/)] ^y $( a c() 
V K n (y,y) J 3 K J V ; ' 

where the O notation is uniform in x, y and u := y — x . 

Proof. We use Lemma 6.1 of [H] which states that if x G (a + e , b — e ), then 

V , f A sin7r(£-?7) 



1 



-,X + 



n\E'(x) " y ' n\l/(x)' ' ri^ix) J ir(^ — rj) 
where the O notation is uniform in x, £ and rj. Thus we have 



fT n (x,x) K n (x,y) 
K n (y,x) K n (y,y) 



K n (x, x)K n (y, y) - K n (x, y)K n (y, x) 
2 sin 2 (7mu\E'(a;)) 



-K 2 n 2 u 2 ^ (x) 2 



0{n) 



-Ti 2 m(xYn\ 2 + 0(n). 



This completes the proof. 



□ 



Lemma 5.2 Uniformly for x, y G (a + eo, b — eo), K n (x, y) = 0(n). Moreover, the first 
and second partial derivatives of K n (x,y) are of order 0(n 2 ) and 0(n 3 ), respectively. 
Finally, if in addition \x — y\ > 5 for some constant 5, then K n (x,y) = 0(1). 

For the proof of Lemma [5.21 see Lemma 10.1 in [I]. 

In conclusion, we successfully extended the result to the UUE case. 



6 Conclusion 

The minimum gap of the eigenvalues is the fine structure of the spectrum. This paper 
addresses two questions: (1) under what scale can we observe a clear picture of them? 
and (2) under the correct scale, what can we see? 

For the Ginibre ensemble where the eigenvalues lie on the two dimensional complex 
plane, the correct scale of the minimum gap is (9(n~ 3//4 ). On the other hand for the 
Wishart ensemble and the universal unitary ensemble where we only have real eigen- 
values, the correct scale is 0(n~ 4 ^ 3 ). We also showed that all of the three cases have 
Poissonian limit after the correct scaling. This implies, heuristically, that the small 
eigenvalue spacing exhibit some asymptotic independency. 

Compared to Vinson's result in j3] and Soshnikov's result in [5], we can obtain the 
joint distribution of k smallest gaps, without requiring k — 1. Our further research may 
include the generalization of the theory into general random point fields, or even in the 
high dimensional case. 
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